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1 Introduction 



In recent years, the study of field theories in a thermal environment has become increasingly popular 
1^]. The interest in this field has increased for a number of reasons. From the experimental side, the 
study of strong interactions at high temperatures and densities has been taken up by means of heavy-ion 
collisions. In order to understand theoretically the outcome of such experiments, we have to learn how 
to treat field theory in a hot and dense environment. Thermal field theory also plays an important 
role in astro-particle physics. Questions concerning the physics of the very early universe can usually 
only be addressed if the impact of the temperatures on the properties of matter is taken into account. 
A famous and in recent years extensively discussed question concerns the phase-transition that might 
have been associated with the restoration of the electroweak symmetry]^ at temperatures of the order 
of 100 GeV. Another interesting phase-transition is the transition from the quark-gluon plasma phase 
of quantum-chromodynamics to the low energy phase where confinement applies and the (approximate) 
chiral symmetry of the standard model is broken. This transition has taken place at temperatures of 
the order of 100 MeV and is expected to be experimentally accessible with the next generation of heavy- 
ion colliders. It is also an interesting open question of particle physics applied to astronomy whether a 
quark-gluon plasma is realized in the interior of neutron stars. 

A large number of open problems is connected to non-equilibrium phenomena. In the physics of the 
early universe this concerns for example the question of reheating after a (hypothetical) inflationary phase. 
Also heavy-ion collisions may not yield a thermalized state and one should strictly speaking consider QCD 
out of thermal equilibrium. Finally it is of interest to study the dynamics of phase-transitions in various 
contexts. 

In marked contrast to the number of situations in which thermal (or more general non-equilibrium) 
field theory is of relevance is the number of open fundamental questions. This is due to the fact that 
in many situations even in the simpler equilibrium case the well established methods of perturbation 
theory fail due to a modified infrared-behavior of the theory. Famous examples of the infrared problems 
are the observation of Linde Q, that the free energy of a non-abelian gauge-theory is not computable 
in perturbation theory beyond three-loop order (irrespective of the size of the coupling constant) or the 
fact that even super-daisy resummed one-loop perturbation theory fails to predict the correct critical 
behavior of self-interacting scalar theories even if they are weakly coupled at vanishing temperature 
[Q]. Although the physical mechanisms behind the above-mentioned problems are quite different, both 
problems may be traced back to the fact that the behavior of bulk-observables like the free energy 
or of theories close to a second order phase-transition is governed by three-dimensional physics (if the 
underlying zero-temperature theory is formulated in 3-1-1 dimensions). 

We then have to devise methods that are able to cope with this effective change of the relevant 

degrees of freedom. One such method in principle is lattice field theory, which may be used to study 

time-independent quantities like the free energy also at non-vanishing temperature. Lattice studies of field 

theory at high temperature are however often complicated by the existence of widely separated scales and 

there is as yet no formulation of lattice field theory that can deal with finite density or time-dependent 

^We know by now that in the framework of the standard model there was no phase-transition [^]. This does however not 
exclude phase-transitions in extended models like the minimal supersymmetric standard model. 
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quantities. 

Another way of dealing with theories where the relevant degrees of freedom are scale-dependent is 
provided by the renormalization group (RG). This is a well known fact in statistical physics, where the 
Wilsonian form of the renormalization group is widely used in the study of critical phenomena. In the 
context of field theory, there are different implementations of the general RG idea. These range from the 
"environmentally friendly" RG (Q and references therein) which is a variant of the Callen-Szymanzik 
RG, via the "auxiliary mass method" to continuum implementations of the Wilsonian approach, ini- 
tiated in a field-theory context by Polchinski and others [^. This approach has in recent years been 
applied to a number of problems in field theory both at vanishing and non-vanishing temperature and is 
now also known under the name "Exact renormalization group (ERG)"-approach. The ERG-approach 
is a nonperturbative formulation of field theory and by construction avoids any infrared-problems ir- 
respective of the dimensionality of the system under study. It is formulated in Euclidean space and 
may be straightforwardly adapted to study thermal field theory in the "imaginary-time"- or Matsubara 
formalism. 

Recently, an implementation of the Wilsonian renormalization group in the real-time formulation of 
thermal field theory has been proposed |§ . This "thermal renormalization group" (TRG) has a number of 
advantages as compared to the ERG applied to thermal field theory. Even though the formulations of field 
theory in the imaginary- and the real-time formalism are in principle equivalent and one can compute any 
quantity of interest in both approaches, in practice recovering Green- functions at real time arguments 
from their expressions at discrete imaginary energies through analytical continuation is often tedious. 
Moreover, in practical applications of the Wilsonian RG results will usually be obtained numerically, 
making analytical continuation impossible. Thus, in order to study quantities like damping rates and the 
like Ijl^, a real-time formulation of thermal field theory is mandatory. 

Also, in the real-time formalism there is a clear separation between "thermal" and "quantum" fluctua- 
tions. The thermal renormalization group only treats the thermal fluctuations of the theory, all quantum 
fluctuations are assumed to be integrated out, i.e. one starts with the full physical effective action of the 
theory at vanishing temperature. This feature essentially prevents us from studying theories that are 
strongly interacting already at T = and may thus seem like a major shortcoming. On the other hand, in 
situations where the zero-temperature theory is weakly coupled this apparent shortcoming turns into an 
advantage. In such a situation the calculation of non-universal quantities such as critical temperatures 
in terms of the measurable couplings of the theory in the framework of the ERG requires a two step 
procedure, where one first has to relate physical quantities at zero temperature to the parameters of the 
action. These parameters are renormalization scheme dependent and one has to perform perturbative 
calculations using the ERG- formulation In situations where two- loop accuracy is required to fix 

the ambiguities of such a calculation, this procedure is highly nontrivial and has up to now only been 



performed in simple models [12|. 

Another advantage of a formulation of RG-equations only for the thermal modes is the fact that 
this can be done while respecting manifest gauge-invariance [13|. The general formulation of Wilsonian 
RG-equations relies on a separation of "hard" and "soft" modes by means of an external scale which acts 
as a momentum cutoff. This procedure is in general not consistent with gauge-invariancepl (see however 



^ There is a formulation using the background- field approach where one can keep manifest gauge-invariance with respect 



3 



1 15 1 for a different approach circumventing this problem), and one obtains generalized Slavnov- Taylor 
identities that may be interpreted as fine-tuning conditions that have to be respected in order to regain 



BRS-invariance of the physical effective action |1(:]. In the framework of the TRG one avoids this problem 
since only the physical fields have thermal fluctuations. Manifest gauge-invariance severely restricts the 
form of possible contributions to the effective action and thus helps in making sensible approximations 
(see the discussion below). 

In the present paper we present an extensive application of the TRG to the critical behavior of 
the simplest nontrivial field-theories in 3 -|- 1 dimensions, scalar self-interacting models with a (at zero 
temperature spontaneously broken) 0(A^)-symmetry. Despite their simplicity, already these theories are 
not accessible in straightforward perturbation theory for temperatures close to the critical one. 0{N)- 
symmetric scalar theories have a number of applications both in the context of statistical physics and in 
particle physics. In statistical physics, these models are used for example for the description of polymers 
(in the limit = 0), the liquid- vapor transition (N = 1, the Ising model), the transition of helium to a 
superfluid state {N = 2) or ferromagnetic systems (the Heisenberg model, N = 3). In particle physics, 
the model with A = 4 at zero temperature describes the scalar part of the Lagrangian of the standard 
model. At non-vanishing temperature the theory with one scalar field shows the same universal behavior 
as the standard model for the critical value of the zero-temperature Higgs mass |17]. For N = 4 the 
universal behavior should be the same as that of the chiral phase-transition in the Nj = 2 chiral limit of 
QCD. 

The universal behavior of the theory may be studied directly in three dimensions. There are accurate 
results on critical exponents from methods like (Borel-resummed) e-expansion, from perturbative calcu- 
lation at fixed dimension, from resummed high-temperature series or from lattice studies (for a review, 
see ||T8|| ). There are also a number of studies using the ERG directly in three dimensions (for example 
Ip^ , pO| , |2l|] ). Also the scaling equation of state (EOS) has been widely studied in the three-dimensional 
theory ([|^]-|2|], see also 0). 

On the other hand, the number of studies of the full 3 + 1-dimensional theory at finite temperature 
close to the phase-transition is rather limited. Results on the critical temperature and exponents in 
the ERG-approach are given in |25]-[28|, and for = 1 Tc, exponents and critical amplitude ratios 
have been studied by means of the environmentally friendly RG in |^9[. The critical temperature and 
some exponents have also been obtained through the auxiliary mass method in |^0|. Results for critical 
exponents from the TRG have been given for A = 1 in and critical exponents as well as the equation 
of state for A^ = 1 have been obtained from the TRG in [ 31 1 . We here extend these studies to arbitrary 
A and furthermore give critical temperatures, the scaling form of the equation of state as well as critical 
amplitudes in terms of the zero-temperature couplings for the cases A = 1 and A = 4. For reasons to be 
discussed below our results on the universal critical behavior are not as accurate as the existing results 
obtained in the effective three-dimensional theory. On the other hand, the present work does not rely on 
universality arguments and the approximations made in the present paper may be improved. The TRG 
allows for the direct computation also of real-time quantities and offers a convenient nonperturbative 
tool to study field theory in a hot and - after a straightforward extension of the method possibly dense 
- environment. 

to background-gauge transformations 
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The format of the paper is as follows: In the next section we give a review of the formulation of 
the thermal renormalization group-equation and discuss some general points regarding the possibilities 
of approximately solving equations of this general type. Here we can rely on experience gained from the 
study of the Euclidean exact renormalization group-equation. In section 3 we apply the method to 0{N)- 
symmetric scalar theories. After deriving the flow-equation for the effective potential we discuss in detail 
the phase-transition for = 1 and N = 4. We compare our results for the critical temperatures, critical 
exponents and amplitudes and the scaling equation of state with different results given in the literature. 
We also discuss the dependence of the critical temperature on the number of fields and compare with the 
expectations from large-A and naive perturbation theory. 

Section 4 contains our concluding remarks. 

2 The Thermal Renormahzation Group-equation 

In this section we discuss the derivation of the thermal renormalization group (TRG)-equation and 
some more general points on its relevance and strategies for its approximate solution. 

The basic objective of thermal renormalization group-equations in the framework proposed in 
is to use the Wilsonian RG-approach in the real-time formulation of thermal field theory in order to 
have access to real-time correlation functions avoiding the need for analytical continuation. Wilsonian 
renormalization group-equations are generally not very useful in Minkowski-space, since the connection 
between k'^ and "softness" of modes is lost. On the other hand, in many situations the problematic 
behavior of perturbation theory at finite temperature is only due to the thermal modes, and the dynamics 
of the theory in the vicinity of a phase-transition is governed by three-dimensional classical statistics. 
To deal with this particular problem, one may use the fact that in real-time the propagators clearly 
discriminate between "thermal" and "quantum" fluctuations to treat the thermal contributions with 
renormalization group methods, leaving the quantum contributions untouched. 

2.1 Formulation of the TRG-equation 

Field-theory, not necessarily in a thermal equilibrium situation, may be formulated using the Schwinger- 
Keldysh or closed time-path (CTP) formulation |^^. Here one is interested in expectation- values of 
operators with respect to some general density-matrix rather than in scattering amplitudes, i.e. matrix- 
elements of the form (out | in) as in the more conventional approach to field-theory at vanishing temper- 
ature. In the CTP-formalism one thus considers matrix elements of the form (in|in). In order to be able 
to derive correlation- functions at different times, one defines the generating functional Z which should 
now depend on two sources (in the case of one real field) Ji and J2 by inserting a complete set of states 
at some time t as 

Z[Ji,J2] = j 5«>(in,to|^,t)j2(«',t|m,to)ji (1) 

where all states are in the Heisenberg picture and Z of course depends on the initial state of the sys- 
tem. One may as usually introduce the path-integral representation for the transition matrix elements 
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Figure 1: The Schwinger-Keldysh contour in the complex time-plane. The lower 
part of the contour is infinitesimally shifted below the real axis. 



{m,to\^,t)J^ to find 



Z[Ji,J2] = J 5^{m,to\T*e^p(^-i J^ dtd^xJ2(l)H^\<^,t) 



x{^,t\Texp ( i / did xJi(l)H ) |in,to) 

'to 



(2) 



where T and T* represent the time- and the anti-time-ordering operator respectively and (j)}{ are Heisen- 
berg-fields. Note that the sources Ji and J2 have to be different in order to obtain time-dependent 
correlation functions. 

Finally, one may also formulate the generating functional by introducing the density matrix p that 
characterizes the initial state. The generating functional then simply is the ensemble-average of the 
product of the time- and anti-time-ordered exponentials and we can write 



Z[Ji,J2,/9] = NTrlp 



r*exp 



-I I dtd xJ24>H 

'to 

t 



-I 



Texp ( i I dtd^xJi(pH 
to 



J to 



dtd-^xJ2(pH \^,t) X 



/ 

J to 



x($,i|rexp i / dtd'^xJi4>H,to {(I), to) 



(3) 



Derivatives of Z with respect to the sources thus generate ensemble-averages (expectation values) of 
products of Heisenberg-fields. Due to the introduction of the two independent sources and the fact that 
we are computing matrix-elements between "in" -states at equal time, while inserting a complete set of 
states at a different time t and using time- and anti-timc-ordcrcd exponentials the integration-contour is 
a closed path from to to t and back, thus giving rise to the name "closcd-timc-path" formulation (figure 
1). Finally we go over to a path-integral representation of the transition elements and write 

z[Ji,j2,p] = j S4>iSMM\p\h,to)e^p[i{s[(t>i]-s*[(j)2] + Mi-J2h)] (4) 

Up to now, the initial state is completely unspecified. Indeed, the matrix element (i?!'i|p|^2) will be 
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some functional of the configurations (piix^t = tg) which may be written as 



{(t)lito\p\4>2,to) = exp (iK [(/)!, (/)2]) 



exp 



(5) 



where we use a compact notation with a,b = 1,2 and metric Cab = diag(+l,— 1). We may then treat 
the coefficients K""^ " as sources for the corresponding composite operators and write the generating 
functional as 

1 



ZIJ^^K, K'',K''\...]=j5(j)i5<j)2exp 



(6) 



Having specified the density matrix, i.e. the sources K""^"', one may then proceed along the lines of 
in order to obtain e.g. a perturbative series etc. 

We will in the present paper be interested in systems with a density matrix of the general form 

d^k 



Pf = C exp 



-f{\k\)alak 



(7) 



(27r)32u;fc 

where and a are creation- and annihilation-operators of a free scalar field respectively. (0) is a 
slight generalization of the corresponding thermal density-matrix which obtains for = PuJk (we 

may interpret as the density-matrix of a system in which each mode of the field is at a different 
temperature given by = — we will come back to this interpretation below). For the above 

density-matrix the matrix-elements of {cl)i\p\(j)2) are given by [p4| ] 



d^k 



[4>l + 4>l) cosh/(|^|)-20i(/<2 



(8) 



(2vr)3 sinh/(|fc| 

We observe that as well as K°'^'^ and all higher derivatives of K[(j)\ vanish in this case. The constant 
is absorbed in the normalization of the generating functional. The fact that only terms quadratic in (pa 
appear in the density-matrix element is of course the reason why in finite temperature field-theory in 
equilibrium only the two-point functions are modified. Indeed, noting that K"'^ only has support at t = Iq, 
we may absorb the effects of the initial state into the boundary-conditions for the free field-equations, 
where they simply modify the two-point functions, yielding in the case of a thermal initial state the 
well-known real-time propagators of the CTP-formalism at finite temperature. 

Above we have made use of the path-integral approach to the CTP-formalism. In simple theories, a 
convenient way to obtain explicit expressions for the propagators is through the operator-formulation. 
To this end, introduce a generalization of the 0-function with real time-arguments to a 0-function on the 
contour C corresponding to eq. dH) and evaluate the path-ordered product (Tc here is the path-ordering 
operator) 



Ocit - t'){4>{x)4>{x')) f + Ocit' - t){4>{x')4>{x)) f 



after substituting for the field operators 



d^k 



-ikx 



+ aie 



(9) 



(10) 
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In d), the averages are ensemble-averages with respect to the density-matrix defined in (0) and in the 
usual manner one obtains, using 

1 



- 1 



5{k — p) 



1 



5{k-p) 



(note that these expressions again reduce to the well-known results for /(|A;|) = f3L0k) 



(11) 



Gf{x-x') 



iDf{x-x') 



-27r(5(A;^ - m^)e 



-ik{x—x') 



Q{ko)Q,{t - t') + @{-ko)Q,{t' -t) + ( exp(/(|fe|) - 1 



(12) 



Returning to the notation in terms of two fields and (f)2 and noting that the time-contour starts at 
Iq ^ —oo and goes to t ^ -|-oo and back shifted by an infinitesimal amount below the real axis (see 
figure 1), and thus ©c obeys 

0c(t — t') = Q{t — t') on the upper part of C 

= Q{t' — t) on the lower part of C 

= for t on the upper and t' on the lower part of C 

= 1 for t on the lower and t' on the upper part of C (13) 

we can write the propagator as a 2 x 2-matrix in the form (note that our sign-conventions are slightly 
different from those in Bl) 



Df{k) 



Ao + (Ao- AS)iVj(|A:|) 
-(Ao-A5) (Q{k^)+Nf{\k\) 



■(Ao-AS)(^e(-A:o)+iV/(|A:| 
-AS + (Ao-A5)iV;(|fc|) 



(14) 



where we have already switched to momentum-space and use 

1 



and 



Ao(A;' 



Nf{\k\) 



A;2 



1 



ojm _ 1 



(15) 



(16) 



We are now in the position to formulate the renormalization group-equations. As stated above, the 
general idea of Wilsonian or exact renormalization group (RG)-equations is to introduce an external 
dimensionful parameter and divide the path-integral into a part that includes only the hard modes and 
a remainder including the soft modes, where the separation is done with respect to the external scale. 
The result of the path-integral over the hard modes is then treated as an effective action. Instead of 
performing any of these path-integrals, one derives functional differential equations for the dependence 
of the effective action T\ on the external scale (we shall denote it by A) and solves these equations for 
A — > 0. In this way one recovers the full solution to the path-integral and thus the solution of the theory. 
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Since at any point one integrates only over an infinitesimal "shell" of momenta between A and A + dA, 
the resulting equations are formally one-loop. 

We have also pointed out that the formulation of exact renormalization group-equations for a theory 
in Minkowski-space at vanishing temperature is possible but not very useful. The reason is that the 
Lorentz-invariant separation of modes with k'^ > (<)A^ as being hard (soft) has no well-defined meaning. 
This translates to the fact that we can give no simple starting- value for the solution of the RG-equations 
as A — > oo. In contrast, in Euclidean space and for a conveniently defined A-dependent effective action 
the boundary condition is simply that Fa becomes the renormalized classical action of the theory as 
A — > oo |35|. To treat field theories in Minkowski-space with Wilsonian RG-equations, we would have to 



use cutoffs that break Lorentz-invariance. 

The situation is different at finite temperature where Lorentz-invariance is broken anyway by the 
presence of the heath-bath with a preferred frame of reference. We could thus introduce cutoffs that, 
e.g., treat ko and k differently. In an imaginary-time formulation this strategy is pursued in ||2^. On the 
other hand, most of the problems of perturbative calculations in thermal field theory are actually related 
only to the infrared-behavior of the thermal modes. Since there is a clear separation in the two-point 
function between thermal and quantum contributions (cf. eq. (|l4|)) and the thermal fluctuations are on- 
shell (note that as e — > we have Aq — Ag —2'7Ti6{k'^ — m'^)), and thus are effectively depending on the 
(Euclidean) space-components of k alone, we may treat these fluctuations with renormalization-group 
methods. 

We proceed as follows: As noted above, with the general form of the density matrix given in (0) 
we may control the boundary conditions for the fields. We have seen that thermal boundary conditions 
correspond to 

f{\k\)=Pu;k = P\/\k\^ + m^ (17) 

If we want the "hard thermal" modes with \k\ S> A to be in thermal equilibrium, we should thus introduce 
the external scale into the function / such that 

fj,{\k\) ^ (3uJk as|fc|»A (18) 

On the other hand, we note that the thermal contributions in the two-point function ( [T^ ) are absent in 
the limit /(|A;|) oo (in the thermal case this would simply correspond to T — > 0). Hence we suppress 
thermal fluctuations with ^ A by requiring 

fA{\k\)^oo as|fc|<A (19) 

Having thus introduced the external scale A we now proceed to derive the renormalization group- 
equations. The easiest way to do this is to turn back to a path-integral representation of the generating 
functional in the presence of the modified density-matrix. Using the inverse of the propagator-matrix 
Df, which we now denote as D^^ to stress the dependence on the external scale. 



j^-i^^ AS-(Ao-AS)iV/(|A;|) -(Ao-A*o)[e{-ko) + Nf{\k\)^ 

^ A0A5 1^ _(Ao-A5)(0(A;o) + A^/(|fc|)) -Aq - (Aq - AS)iV/(|fc|) ' 
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we have 



)id(p2 exp 



i ( (^A^)''^ h + Sint[01, 02] + J>a 



(21) 



where S'int is the interaction-part of the classical action with S\at[4>i-,4>2\ = 5'int — 'S'int[02]- It is 
straightforward to derive an equation for the dependence of Za on A. One finds 0] 



-'-Tr 
2 



(22) 



where here and below we use the compact notation J = (Ji, J2) and correspondingly for the fields. 
Defining as usually the generating functional for the connected Green functions Wa[J] = — ilnZA[J] we 
obtain the RG-equation for W\, 



-Tr 



d 



+ -Tr 
2 



'5Wa[J] f^±jj^i\ ^Wa[J] 



6J 



6J 



5J5J 

Finally, we define the generating functional of IPI Green functions as the Legendre transform of W\: 

1 / „ _ X \ ab 



(23) 



where 



TaM = Wa[J] - - ^Va [Dl^) 

5VFa[J] 



5J« 



(24) 



(25) 



and we have subtracted a term bilinear in the field, corresponding to the free (cutoff-)propagator. In the 
RG-equation for Fa, subtracting the free propagator leads to a term that cancels the contribution from 



the second term on the right hand side of ( 23 ) and we obtain the thermal renormalization group-equation 



(26) 



Finally we have to discuss the limits A ^ cxd and A — > 0. Let us start with the limit A ^ 0. In this 
case, the density-matrix according to (0) and ( [l8|) is the usual density-matrix for a system at thermal 
equilibrium with temperature (3^^ ■ Accordingly, the inverse propagator becomes the free thermal 
propagator (cf. eq. (pO|)) and the effective action tends to 



|irn FaM =T[ip\--ipa {D 



-1 



, ab 



(27) 



where ( |24D was used. 

On the other hand, to consider the limit A 00, note that limA-+oo = according to (|T^ ) and 
(p!9|). In this limit, the propagator is the one of the theory at vanishing temperature (cf. eq. (|14|)). 
The starting value of the flow-equation as A ^ 00 is thus the effective action of the CTP-formalism at 
vanishing temperature. 

Before proceeding to an application of the formalism, let us add some general remarks. First, note that 
even though equation (|2^) is an exact result, it may be formally considered as a one loop-equation. The 



10 



^^~9A~ - 2 



Figure 2: Graphical representation of the thermal renormahzation group-equation 
&■ 

corresponding diagram is given in figure 2, where the squared cross denotes the insertion of Ad/^D]^^ and 
the hne with the dark dot is the full propagator. Of course, in an interacting theory the full propagator is 
still a functional of the fields, making the equation impossible to solve exactly - nevertheless the existence 
of a pictorial description in terms of Feynman-like diagrams simplifies the discussion of the flow-equations 
and allows us to easily derive flow-equations for Green-functions etc. 

One should note that in the context of thermal renormahzation group-equations as in the usual real- 
time formalism one cannot generally dispose the thermal ghost fields (992 in the above derivation). Even 
if one is only interested in Green- functions with physical (cpi) external lines, the facts that equation 
(|26| ) is formally of one-loop order and that the classical vertices do not mix the fields do not help since 
taking derivatives of (|2^ ) with respect to tpi introduces full vertices with external legs. These of course 
will in general mix physical and ghost-fields. Also note that AOaDJ^^ is a non-diagonal matrix. Since 
the insertion is basically treated like a vertex that mixes physical and ghost fields in ( p^ one gets 
contributions from all matrix-elements of the full propagator. 

Finally, from the derivation given above it is immediately obvious that the general method of thermal 
renormahzation group-equations is not restricted to thermal field theory in equilibrium. In principle one 
may use any form of density-matrices for the initial state. It is clear, however, that if the density-matrix 
is not bilinear in creation- and annihilation operators there will be A-dependent contributions also to 
the interaction terms in the path-integral (^l]) and the renormahzation group-equation would not have a 
simple form. It is also straightforward to formulate thermal renormahzation group-equations for gauge 
theories | ]13| or theories involving fermions. 

Apart from the fact that the formulation of renormahzation group-equations in the real-time approach 
to thermal field theories simplifies the calculation of real-time correlation functions since one avoids ana- 
lytical continuation of quantities calculated at imaginary time arguments, the major difference between 
the thermal renormahzation group formalism reviewed here and the Matsubara-approach by means of 
exact renormahzation group-equations is the boundary condition for A ^ 00. As we have seen, the 
starting value for the effective action in the TRG-formulation is the full effective action at vanishing 
temperature. This should be contrasted to the imaginary-time approach where one treats quantum- 
and thermal fiuctuations on equal footing. There the starting value for the effective action is the renor- 
malized classical action. In practice this means that for a computation of finite-temperature quantities 
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in terms of physical quantities at vanishing temperature a two-step procedure is necessary pq ]. One 
first has to connect the parameters of the classical action to observables in a T = 0-computation^ and 
only afterwards calculate quantities at finite temperature. To get control over the relation between e.g. 
perturbatively defined MS-parameters and the corresponding couplings in the renormalization-scheme 
defined by the exact renormalization group approach one has to perform perturbative calculations also 
in the RG-approach In cases where one needs to go beyond one-loop order, this is an extremely 

cumbersome task jl^. In this respect the approach used in the present work is very convenient, since 
perturbative contributions at vanishing temperature may be taken into account by using perturbative 
starting conditions to any desired order. This feature implies of course an obvious limitation to the use 
of thermal renormalization group-equations in situations where already the theory at T = is strongly 
coupled. In such a situation one needs input from other nonperturbative methods such as lattice com- 
putations or exact renormalization group-equations at vanishing temperature. On the other hand this 
is a very welcome feature if one wants to study questions like the validity of perturbative dimensional 
reduction, a method that e.g. is at the heart of most of the present nonperturbative results on the 
electroweak phase-transition in the early universe Q and is also studied in the context of QCD (for a 
review see e.g. fs^). In the case of the electroweak phase-transition there are a number of lattice studies 
that do not rely on perturbative dimensional reduction at least in the bosonic sector of the model (see for 
example and references therein) . A comparison of these results with studies assuming the validity of 



perturbative dimensional reduction indicates that method works fairly well in this sector [38|. However, 
to date there exist no nonperturbative checks of the method in the fermionic sector of the theory which 
is quantitatively important due to the large top- Yukawa coupling. TRG-equations should allow a test of 
this method also for models with chiral fermions. 

2.2 Approximation schemes 

Before presenting the results of our study of self-interacting scalar theories in the next section, we would 
finally like to briefly discuss some general strategies to approximately solve equations such as ( p6D |1S]. 



The full effective action T\ may be characterized by infinitely many couplings multiplying the invariants 
consistent with the symmetries of the theory under study. The functional differential equation (^6|) is thus 
equivalent to an infinite system of coupled nonlinear ordinary differential equations. Clearly one has to 
make some approximation to proceed. All approximations involve a truncation of the system of ordinary 
differential equations to another (possibly still infinite) system of equations by neglecting couplings. 

The simplest approximations reduce the number of couplings considered to a finite set. A prominent 
example of such a truncation is the so-called "local polynomial approximation" (LPA) to the effective 
potential, supplemented by a derivative expansion. This approximation is at the heart of most applica- 
tions of Wilsonian RG-equations considered in the literature. It may be viewed as an expansion in the 
canonical dimension of the couplings (see [^]). 

Alternatively one may consider truncations that keep an infinite set of couplings. This will result 

in (systems of) partial (integro-)differential equations. There are two more or less orthogonal ways to 

^The exact renormalization group-equations imply a specific renormalization scheme. One may thus not directly compare 
results from perturbative computations in, say, the MS-scheme with renormalization group results at the same values of the 
parameters of the renormalized action. 
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proceed: 



(i) One may perform a derivative expansion of the effective action. In this case one classifies the 
invariants by the number of derivatives appearing. The coefficients of these invariants are then 
functions of constant fields. Considering a scalar theory with one field (p and Z2-symmetry for 
simplicity, one expands the effective action as 



(fx 



(28) 



where we have taken into account the fact that at finite temperature Lorentz-invariance is broken by 
the rest- frame of the heat-bath and have introduced the corresponding four- vector u^. Truncating 
the effective action at finite order in this expansion, one derives a set of coupled partial differential 
equations for the functions VA^cp^), Zj^^tp^) etc. 

This is the approach that we will pursue below. 



(ii) One may expand the effective action in powers of the fields, 

1 /" ,^ in) , 



rA = Ta 



^=0 + X] 

n=l 



—.Id Xi...d XnTh [Xl, ...,Xn)'p{xi)...(p{Xn) 



(29) 



The coefficients of this expansion are the IPX Green-functions (the full propagators and vertices) of 
the theory. This is the strategy that is naturally followed in the study of Schwinger-Dyson equations 
13£ 



The approach using the derivative expansion is technically relatively straightforward and has been 

21] and at finite temperature [31, matrix- 



used to study scalar theories in three dimensions 



models |41], the abelian as well as an S'C/(2)-Higgs model in three dimensions |42|, and an effective 
theory for low-energy QCD at finite temperature in the imaginary-time formalism [^]. It is generally 
appropriate in situations where the degrees of freedom considered correspond also to asymptotic states 
and the anomalous dimensions are small. As noted above one obtains a set of coupled partial differential 
equations for functions that depend on the invariants that one may construct from the fields considered 
(at vanishing momenta). 

On the other hand, an expansion in the n-point functions rapidly becomes technically involved with 
growing n. Considering only the flow-equation for the two-point function of a scalar theory at T = 0, 
one has to deal with a partial integro-differential equation of the structure 



d 



+ 



dH 
d'^l 



Fa 



+ 



(27r)' 



(2)^;2^ .,(3). 2 ,2 



(30) 



corresponding to the diagrams depicted in figure 3. In a theory with vector-fields like QCD the number 
of invariants that make up a specific n-point function is rapidly growing with n. So far, this approach 



has been used to study the 2-point functions of pure Yang-Mills theory [44| and the four-point function 
of an effective model for QCD pH] (both at vanishing temperature). 
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P 

Figure 3: Exact flow-equation for tlie two-point function in a scalar theory. 

Which of the above mentioned approaches should be considered depends on the physics of the theory 
under study. Generally, the derivative expansion will not be useful in situations where the momentum- 
dependence of the correlation functions is strongly modified by interactions. 

After this more general review we will in the following section consider self-interacting scalar theories 
with an 0(A^)-symmetry at finite temperature. 




3 Critical behavior of 0(A^)-symmetric scalar field theory 

0(A^)-symmetric scalar theories are well studied both in 4 dimensions at vanishing temperature as well 
as in statistical physics. The 0(A^)-symmetry of the action may be spontaneously broken at zero- 
temperature. In this case, one expects symmetry-restoration for large temperatures. The universal 
aspects of the corresponding phase-transition are very well known from considerations of the model in 
three Euclidean dimensions (for a review, see ||l^). Quantities like the critical exponents are calculated 
to high orders in the e-expansion or in perturbation theory at fixed dimension (i = 3, where in both cases 
Borel-resummation has to be used to deal with the divergent series. There are also results from lattice 
studies available. 

What is less well studied are the non-universal aspects like the critical temperature as a function of 
the zero-temperature "physical" couplings. To study non-universal quantities one really has to consider 
the model at finite temperature. Perturbative studies of the theory at finite temperature are complicated 
by the bad infrared behavior near the phase-transition. Due to this behavior, for = 1 even super-daisy 
resummed perturbation theory falsely predicts a first-order phase-transition Q. Alternatively, one may 
employ an expansion in 1/A. In this case the resummation of daisy-diagrams yields exact results in 
leading order in 1/A [^] (see also |l^). One finds a second order phase-transition with mean-field values 
of the critical exponents /? and r/, whereas v = 1 + ©(l/A), 5 = b + 0{l/N), a = -I + 0{l/N) and 
7 = 2 -|- 0{\/N) (see below for definitions of the critical exponents). 

As noted in the introduction, there are a number of methods that correctly predict a second-order 
phase-transition of the finite-temperature theory also for A^ = 1. Results for the critical temperature 



have for example been given in [g^, 29] and we will compare the corresponding values below. 



3.1 Flow-equation for the effective potential 

As stated in the last section, there is no way to solve the TRG-equation ( ]26| ) exactly and we have to make 
simplifications. In the theory under consideration here, many studies have shown that a good truncation 
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consists in performing a derivative expansion of the effective action and considering only the first few 
contributions. This fact is related to the smallness of the anomalous dimensions of the fields, which 
govern the scale-dependence of the wave- function renormalizations. In this spirit we will in the following 
approximate the effective action Fa by its value for constant configurations, the effective potential Va, 
and use a standard kinetic term. The effective potential is of course still a function of the fields. We 



may at this point use arguments given by |47| to simplify the task of calculating the effective potential 
considerably. The line of argument goes as follows: The free energy of the theory is given by the functional 



T\[ip] which is defined as |47] 



(31) 



6ip 5ipi 

We are interested in the value of the effective action for constant field configurations and define the 
A 



tadpole T^p^{lp) through 



By virtue of the symmetry 



(32) 



(33) 



one may derive relations between n-point functions with type-1 and type-2 external legs. 

In the following we will discuss 0(A^)-symmetric scalar theories. We will label the physical- and 
thermal ghost-fields with subscripts as above, and use greek letters as superscripts to label the components 
of the 0(A^)-vectors {a, (3, ... = 1...N). We frequently use the 0(A'^)-symmetry and write 



Now consider the renormalization group-equation for the tadpole 
respect to cpj and setting the fields constant we obtain 



(34) 

Deriving equation ( p6|) with 



A 



d ST A 

dA 6^1 



(P) 




AoA* 



dA 









;)] 



(35) 



where the trace is over the "thermal" matrix-structure and we have used the fact that the two-point 
functions are diagonal with respect to the 0(A'')-indices. The term in square brackets constitutes the 
"kernel" of the TRG-equation and corresponds to the (modified) logarithmic derivative of the full 
(cutoff-)propagator with respect to A (no sum over a). 



-(Ao 



AqAS 







dA 



(; 






01 



5^Ta 



-1 



d 



(36) 



15 



where di^\ denotes a derivative acting only on the exphcit A-dependence of the free propagator DJ^^ . 

Equation ( p5| ) is an obvious generaHzation of the case of one scalar field studied in 31 ] . The trace 
yields 



tr 



(; 


;)( 






-as: 


AaA* 







(37) 



A 



(a) 



(38) 



where the self-energy is given by (no sum over a) 



Im U^^\k,^) 



Re 



6^f{-k)6^f{k) 



-Im 



V^l=f2 =</3=const 



l + 2Nf{\k\) dipf{-k)6ipf{k) 



(39) 



To the considered order in the derivative expansion, the self-energies are momentum independent and 
real. Then in equation (^) we can replace ^ 



i(Ao-AS) 
and may use 



27r5 lk'-m' + U^^\ip))e[\k\''+m'- U'^^' (if) 



AoA* 

to write (no sum over a) 

n?)(<^) 



(40) 



(41) 



We finally have to specify the form of the cut-off function /a- In order to be able to do the loop-integral 
in ( p5| ) analytically, following 31 1 we will choose this function such that 



(42) 



where A'^^^; denotes the thermal distribution function for bosons. In this case we have d\N\{\k\) 
NBEi^k)^i\k\ — A) and the loop integrals are trivial. 
We write for the functional Fa 



1 



1 



1, 



d^x -mV' - Va{^) + -ZM {d^Y + o^a(9^) (n • d^Y + 



(43) 



and, working to lowest order in the derivative expansion, neglect Z\ and Ya- We will allow for spontaneous 
symmetry breaking along the "l"-direction in the 0(A^)-group, giving rise to a radial (Higgs-) mode, 
the field (/?^, and A^ — 1 Goldstone-modes t/?"^^. This amounts to evaluating the effective action on a 
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background field configuration = i/2p5"^ + x'^il) ^iid putting = at the end. We then have (no 
sum over a) 



ii^ [p)-m + g^^^^^ 



(44) 



Dropping the subscript A we find from ( pjq ) 

d dV{p) 



A 



dA dif-r 



A3 Nbe{ooi) d^Vjp) 2 

4^2 9^79(^15(^1 ^^""1^ 

-(A'-i)^^.ISLe(.i> 



LOi and denote the frequencies of the radial and the Goldstone-modes respectively, i.e. 



(45) 



UJl 



A2 + -^-^ = VA2 + V'ip) + 2pV"ip) ; c^2 = WA2 + 



92 y 



VA2+T%o) (46) 



where primes denote derivatives with respect to p. We may finally integrate eq. (45) with respect to ipH 
to obtain the flow-equation for the effective potential as 

A' 



-T— In 
27r2 



1 - exp [-13^ + V'{p) + 2pV"{p)^^ + V'{p) + 2pV"{p)) 
A3 



-{N - l)T—^ In [l - exp ^-(3^ + V'{p))\ Q{A^ + V'{p)) 



(47) 



We will in the numerical evaluation also use the flow-equation for the minimum of the potential in 
the spontaneously broken phase, po- This is obtained from the condition y'{p)\pQ = and reads 



A 



dpo 
dA 



47r2 



(3)' 



3 + 2 



K ) TA^TWexp JW^^ - 1 



e(A2 + 2poV^) + 



A^ 1 



1 



47r2 A [a? 1 
exp a / ^ - 1 



(48) 



As we will see shortly, it is also convenient to define dimensionless quantities in the following way: 



The corresponding flow-equations read 

1 



u 



V 



dp* dK* 



(49) 



du 



3u + UK K In 

27r^ 



1 - exp -AVl + u' + 2ku" 6(1 + v! + 2kx^' 



iV- 1 



In 



and 



-Ko + 



4^2 



3 + 2 



1 - exp (-AVl + 



(3)' 



9(1 + u') 



(50) 



VI + 2kqvI^ gAVl+2Ko< _ 1 



9(1 + 2koO + 



+(iV-l^ 



A 1 

47r2 - 1 



(51) 
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As was already pointed out in 31], if one considers the limit where A"^ + V' + 2pV" and A'^ + V' are both 
<C one recovers the flow-equations of the purely three-dimensional theory with a sharp cutoff up to a 
field independent contribution. This is the relevant limit in the scaling-regime close to the phase-transition 
and motivates the use of the rescaled variables defined in (|49|). The powers of the temperature appearing 
in eq.s (^) are identical to those obtained by naive dimensional reduction, whereas the powers of A are 
the ones one would use in the purely three-dimensional theory to make the couplings dimensionless. 

We finally have to give the starting value for the effective potential in the limit A — > oo. As discussed 
in section 2, the effective action in this limit becomes the effective action of the zero-temperature theory. 
Correspondingly the starting value for the effective potential should be the full effective potential of the 
theory at vanishing temperature. At this point we face an obvious problem if we are to use perturbative 
results on this quantity. Consider the effective potential of an 0(A^)-symmetric scalar theory in the broken 
phase at T = as computed e.g. in the MS-scheme. Taking the scale to be set by the zero-temperature 
tree-level Higgs-mass uiq = 2gpQ to one loop we have 



Vt=o{p) 



,92, 9 



(3p - pqY 



In 



+(iV-l) 



9 



647r2 



(p - PoY 



3p - Po 

2po 
- jOQ _ 3 
2 



+ 



2po 



(52) 



On the right hand side of the flow-equation ( [471) we need the flrst two derivatives of the potential with 
respect to p at arbitrary values of the field. Now from equation ( |5^ ) it is obvious that the second derivative 
of the potential has logarithmic singularities at p = ^ and, for N > 1, a,t p = pQ. This is of course a 
well known fact which for example forces us to define the renormalized quartic coupling in the massless 
theory with = 1 away from the point p = |48]. As long as we work with a polynomial expansion, i.e. 
as long as we choose to write the potential as 



V{p) 



(53) 



and derive flow-equations for the couplings V^^\ this is no problem. For = 1 we may use p = pQ, 
for A^ > 1 one may still proceed as long as one does not expand the potential about its minimum. 
Even in that case, we could in principle avoid any infrared problems by using the perturbative effective 
potential evaluated with an infrared cut-off in the loop integration. This would compare closely to the 
approach followed in |25], where a local polynomial approximation was used for the potential and the 
zero-temperature renormalized couplings were defined at some non-vanishing scale. In this work we will 
restrict ourselves to small values of the zero-temperature couplings and neglect the quantum-corrections 



to the effective potential, that is we use as a starting value for the integration of the fiow-equation (|47D 
the tree-level potential 



VAoip) 



9 2 

-9P0P + ^P 



(54) 



This simplification has of course no effect on the universal aspects of the critical behavior such as the criti- 
cal exponents, certain ratios of critical amplitudes or the scaling equation of state. For non-universal quan- 
tities such as the critical temperature we will have corrections from the (logarithmic) four-dimensional 
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running of the couplings. These corrections should however be small for the values of the couplings con- 
sidered in the present work. We now discuss the results of a solution of the flow-equations given in (p7[), 
( |48| ) , (|50| ) and (|5^) with the starting value given by (|^) (or the corresponding scaled potential according 
to (^9|)) for different values of N. For the numerical work, the methods proposed in are used. 

3.2 Numerical results 

We start by discussing the results for the case of one scalar field, corresponding to the Ising model of 
statistical mechanics. The first question concerns the order of the phase-transition. As was already 
pointed out above, daisy-resummed perturbation theory fails in predicting the correct critical behavior 
of the theory in this case. In the framework of thermal renormalization group-equations, the case = 1 
has already been addressed in |^, 31], where the transition was found to be second order and critical 



exponents have been given. 

We display the dependence of the minimum of the potential on the external scale A at different 
temperatures in figure 4. In the upper panel, the dimensionless minimum kq is displayed as a function 
of In ^ for several temperatures around the critical one (we have chosen g = 0.1 here). We start the 
evolution of the potential at some large value of A/T in the broken phase, where po{T = 0) is used to 
set the scale and all dimensionful quantities are given in units of po(T = 0). For temperatures larger 
than the critical one (solid lines), the minimum vanishes at some nonzero scale Asymm and the theory 
is in the symmetric phase. On the other hand, for temperatures below Tc the dimensionful minimum 
po approaches some constant as A — > (this is seen in the lower panel of fig. 4, where we display the 
dimensionful minimum as a function of A/T). According to (|49| ) kq then diverges. We see from figure 4 
that at the critical temperature, the dimensionless minimum asymptotically reaches a finite nonvanishing 
value, Kq. As is expected from dimensional reduction and was explicitly demonstrated in [31|, this value 



is the fixed-point value of the corresponding three-dimensional theory. The phase-transition is of second 
order and the theory is in the universality class of the three-dimensional Ising model. Before turning to 
the universal behavior of the theory, let us briefiy discuss the critical temperature as a function of the 
zero-temperature coupling. Our results on this quantity are displayed in figure 5, where we plot the ratio 
of Tc and the naive perturbative result, given by 
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as solid line and compare to the values found in the ERG-approach in the Matsubar a- formalism in |25] 
as squares. We find good agreement as expected^. 

We finally discuss the results on the universal behavior as obtained from the TRG. The universal 
critical behavior of a given field theory is summarized in a number of critical exponents, certain ratios 
of critical amplitudes, and the scaling equation of state (see e.g. [^]). We will start by discussing 
critical exponents and amplitudes. These quantities encode the behavior of different observables as the 
temperature approaches its critical value. Specifically, we will consider the following quantities: 



''Note that the results quoted in were obtained with an exponential cutoff-function rather than with the 0-function 
cutoff used here. 
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Figure 4: Upper panel: The dimensionless minimum as a function of In A/T for var- 
ious values of the temperature around Tp. For temperatures above Tc the minimum 
vanishes at some finite scale Agymm) for T <Tc the dimensionless minimum diverges 
as A ^ 0. Lower panel: The dimensionful minimum po (in units of po{T = 0)) as a 
function of A/T. For temperatures below the critical temperature the dimensionful 
minimum approaches a constant value that is a decreasing function of the temper- 
ature until it vanishes at Tc as A 0. For low temperatures thermal fluctuations 
play effectively no role and po almost remains at the starting value. 
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Figure 5: The ratio Tc/Tc^i-ioop as a function of the zero-temperature coupling g for 



= 1. For comparison results from |25] are also given (squares). 



The (unrenormalized) mass of the order-parameter field (in the language of statistical mechanics 
the inverse magnetic susceptibility in zero field). The susceptibility diverges at the phase-transition, 
corresponding to a behavior of the mass according to 



M| = (C7±) 



Tr-T 



(56) 



This defines the critical exponent 7 and the amplitudes C±. The ratio C^/C^ is a universal 
quantity. 



(ii) 



The quartic coupling g(A 0) goes to zero at the critical temperature according to 

g± = 



Tc 



(57) 



m 



The exponent v also governs the behavior of the renormalized mass at the critical temperature. 
Since we have approximated the wave- function renormalizations to be constant, u should be equal 
to 7/2. 

The dimensionful minimum of the effective potential (corresponding to the spontaneous magneti- 
zation) approaches as the temperature approaches Tc from below, where 

'Tr-T^l^ 



ipo = B 



Tc 



(58) 



(iv) 



Finally, the behavior of the effective potential at the critical temperature is described for small 
values of p by (remember ip = y^2p) 

dV 



dip 



Dip" 



(59) 
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/3 


7 


6 


r? 


V 


This work 
TRG+LPA |] 


0.345 


1.37 


4.97 
3.57 


0.015 


0.67 
0.58 


3d ERG + LPA |1S] 

3d ERG ^] 
3d ERG + WFR ^ 


0.333 
0.336 
0.330 


1.247 
1.258 
1.232 


4.75 


0.045 
0.044 
0.047 


0.638 
0.643 
0.631 


Best values |18| 


0.325 


1.240 


4.81 


0.032 


0.630 



Table 1: Critical exponents from different approaches (A^ = 1) 



We stress that we are considering the theory explicitly at finite temperature. Thus one may expect 
the scaling behavior encoded in the critical exponents only in the region where dimensional reduction is 
effective. The relation (|5^) will for example not hold for large (/? ^ T where the three-dimensional limit 
is not reached. Effective exponents taking into account this "dimensional crossover" are discussed for 
example in |^] and the second reference of |^7|. Our results for the critical exponents defined above are 
given in table 1, the results for some critical amplitude-ratios and critical couplings are displayed in table 2. 
In both tables, we give results obtained by other approaches for comparison. Let us make some comments 
regarding the results presented in the tables. We should first note that the exponents may, apart from 
the possibility discussed above, also be determined from the scaling equation of state to be discussed 
below. We have in particular obtained the value of b in this way. The value is completely consistent 
with (5 = 5 which was found using the thermal renormalization group, but extracting b according to the 
above definition, eq. (|59|) in ||3l|. In fact, one can show that (5 = 5 is the exact solution to the fixed-point 



equations in the three-dimensional theory in lowest order in the derivative expansion |50|. 

The other exponents given in the table have been determined in a number of ways: We have deter- 
mined /3 from ( |58| ) above and from the scaling equation of state and find values of 0.3458 and 0.3451, 
identical within the numerical accuracy. 7 was obtained from the definition given in ( [5^ ) and from the 
asymptotic behavior of the Widom scaling function parametrizing the equation of state and we find 1.369 
from the EOS, 1.368 from (56) if is approached from above and 1.369 if we approach the critical 
temperature from below. Finally v is obtained from (^) to be 0.679 coming from T > Tc and 0.663 
for T ^ Tc from below. All in all the critical exponents are rather robust with respect to details of the 
numerical procedure. The deviation from the best values given also in table 1 is due to the approxima- 
tions made in reducing the exact renormalization group-equation to the flow-equation for the effective 
potential, eq. (47), namely the neglect of higher orders in the derivative expansion. This also prevents us 
from studying the critical exponent rj which is defined as the anomalous dimension at the critical tem- 
perature as A — > 0. From the table it is also clear that we could improve on the results for the exponents 
and amplitudes by taking into account the wave- function renormalizations, as has been done in the work 
listed as ERG in table 1. In p!9| ] the potential was expanded around the minimum and only a finite 
number of couplings was kept (LPA). On the other hand, the authors went beyond leading order in the 
derivative expansion. The LPA was dropped in [pOl, where however the wave-function renormalization 



was still assumed to be field independent. This approximation is finally also given up in [21|. We should 
keep in mind that all these results were obtained in the framework of a three-dimensional effective theory 
rather than in the full temperature-dependent problem. 
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r / 
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'^4* 




B 




D 


(5T=o = 0.1) 


This work 


4.50 


1.76 


19.6 


1725 


1.05 


0.54 


2.73 


3d ERG + LPA |19[ 

3d ERG ^ 
3d ERG + WFR |21] 


4.29 
4.966 


1.61 
1.647 


27.8 


1311 








Monte-Carlo [^| 






23.3 


1476 








Best values |1^ 


4.95 


1.65 













Table 2: Critical amplitude ratios and couplings from different approaches {N = 1). The last three 
entries are non-universal. 



Let us now discuss the results presented in table 2. Here we give the universal ratios of critical 
amplitudes C+ / C_ and as well as the universal critical couplings A4* and Ag* and the non- universal 
amplitudes B, C_|_ and D for gT=o = 0.1. The definition of C± has been given in (^6|). Ry. is a universal 
combination of several non-universal quantities and reads 

R^ = C+DB^'^ (60) 

We have extracted the numbers for the amplitudes used here from the equation of state to be discussed 
below. We have also extracted the ratio C+/C_ from the behavior of the unrenormalized mass as given 
in ([5^). In this case we find C+/C_ = 4.58, consistent with the value quoted in the table. As noted in 
the discussion of the critical exponents, also in the case of the amplitudes the major source of error is 
the derivative expansion used here, as may be seen by again comparing with the results given in |pT|| . 
The universal critical couplings given in the table are defined by 

00 , 

Vi^) = E ^^'^ (61) 

and are taken at the critical temperature, i.e. in the scaling limit. A4* and Ag* are thus simply related to 
the second and third derivative of the potential at the critical temperature with respect to p. The values 
given for comparison have been found through the ERG, using a local polynomial expansion and on the 
lattice, in both cases working in three dimensions. 

Let us then turn to the scaling equation of state. In the language of statistical mechanics the equation 
of state relates the temperature, the order parameter and the external field. In the present setting, the 
external field H is given hy H = V'{p)ip and the equation of state has the following scaling form: 

The function f{x) is known as the Widom scaling function and is universal up to normalization and a 
rescaling of x. It encodes information about the universal critical behavior and several amplitudes and 
critical exponents may be obtained from / in different limits. For the amplitudes and exponents treated 
above we have the following relations involving the equation of state: The value of / at x = gives the 
amplitude D, 

lim f{x) = D (63) 

x—^O 
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Figure 6: The logarithm of the Widom scaling function f{x) for N = 1 (solid lines). 
The left panel shows f{x) in the broken phase, corresponding to x < whereas 
the right panel gives / in the symmetric phase. The short dashed lines are from 
the ERG in three dimensions [20|, the long dashed lines give results from lattice 



simulations (broken phase: [24|, symmetric phase: 



The asymptotic behavior of / in the symmetric phase is 

lim f{x) = {C+y'^x^ (64) 

X— > + oo 

The amplitude B and the critical exponent /3 are related to the zero of / in the broken phase through 

/ (-B-1//3) = (65) 



and the amplitude C- is obtained from the zero of / and the value of the derivative at the zero through 



(66) 



We have pointed out that the normalization of / and of x are non-universal. We will in the following 
figures compare our results with results from exact renormalization group-equations in three dimensions 



and with results from Monte-Carlo simulations in the broken [24| and in the symmetric phase |23] 
(figure 6). In order to do so, we fix the normalization by demanding equality of the results given in form 
of a numerical fit in |20| and approximate polynomial expressions in |2|, with our results for 2 values 
of X. Since the lattice simulations used for the comparison were obtained on different lattices we have 
used different rescalings for the Monte-Carlo results in the left and right panels in figure 6, corresponding 
to the broken and symmetric phase respectively. In order to have a fair comparison also of the results 
from |2C] with the lattice results in both phases, we have also rescaled these results differently in both 
panels. 
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Figure 7: /(x) for = 1 in the symmetric phase compared to the three-dimensional 
ERG results from [pH] for different normalizations (see text for details). 



In principle, for a comparison of our results with the ones of [20| we only have two free parameters 
corresponding to two values of x where equality may be imposed. Thus we also present in figure 7 a 
comparison of our curve with the curve given by |2^] for different normalizations: The solid line in figure 
7 corresponds to our result in the symmetric phase, the short dashed line is the same as in the right 
panel of figure 6, obtained by equating the results at two values of x in the symmetric phase, whereas 



the dashed-dotted line is the result from [20| if we use the two values of x from the broken phase (left 
panel of fig. 6) for normalization for the whole range of x. 

The curves obtained in this work and the results given in ||2^ are almost identical in the broken phase 
(left panel of fig. 6) and show only a slight deviation in the symmetric phase for large x, connected to 
the different results for the critical exponent 7 (see table 1). The agreement with the results from lattice 
simulations is rather satisfactory, as was already noted in |^T| for the broken phase. In view of the sizable 
deviations of results obtained in the framework of the e-expansion or other perturbative approaches from 
the lattice results^, this is a nontrivial achievement. The deviation in the asymptotic behavior for large 
positive X in figure 6 is of no concern. In this region the approximate polynomial expression given in 
p3[ is incompatible with the known behavior corresponding to (64) - it would give an exponent 7 = 2 in 
contrast to the known results (table 1). 

After presenting detailed results for = 1, we now proceed to a study of the case = 4. This case is 
interesting, since the is expected to govern the universal behavior of the two-flavor chiral phase-transition 
in the chiral limit. The arguments in favor of this conjecture rely on dimensional reduction and the fact 
that in the imaginary-time formulation the fermionic degrees of freedom obey anti-periodic boundary 
conditions and thus have no static modes. The infrared behavior should then be dictated by the bosonic 



^For a comparison of the lattice results with results of the e-expansion, see 
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Figure 8: The ratio Tc/Tc^i-\oop as a function of the zero-temperature coupling g for 
= 4. For comparison results from pH] are also given. 



excitations alone, which in this case are the pions and the "sigma" (for a more complete discussion of 
the chiral phase-transition in the two-flavor case also away from the chiral limit in the imaginary-time 
formalism see e.g. [^]). We should also mention that due to the existence of massless degrees of freedom 
also in the broken phase away from the critical temperature perturbative calculations for > 1 are even 
more infrared problematic then for the case discussed above. 

We start our discussion by plotting the critical temperature as a function of the zero-temperature 
coupling g{T = 0) in figure 8 (solid line). We again display for comparison the values obtained by Tetradis 
and Wetterich |2^. As for the case A^ = 1 we find excellent agreement of the results. 

Next we turn to the critical equation of state. This is displayed in figure 9, where as for A^ = 1 we 
display ln(/) as a function of ln(x) in the symmetric phase (x > 0) in the right panel and as a function of 
ln(— x) in the broken phase {x < 0) in the left panel. As discussed above one obtains critical exponents 
and amplitudes from the scaling function, and we collect the values together with some results from other 
approaches in tables 3 and 4. 

Considering the Widom scaling function, we again observe from figure 9 almost perfect agreement of 
our result (the solid curve) with results from the ERG given in [^] (dashed curve) in the broken phase. 
In the symmetric phase, there is a slight deviation for large x, which again is due to the different values 



found for the exponent 7 (table 3). In figure 9, we have normalized the results from [51| to our results for 
two positive values of x close to x = 0. A comparison of the results from the ERG and results from other 
approaches including lattice and e-expansion may be found in |51] and again shows that the lattice results 
on the equation of state compare well with results from the Wilson renormalization group, whereas the 
e-expansion differs from the Monte-Carlo results. 

For the critical exponents given in table 3 we note a slightly larger deviation of our values from those 
obtained from the lattice |52] or from higher order perturbation theory [E^]. Although the values obtained 
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Figure 9: The Widom scaling function for = 4 (solid line) in the broken (left 
panel) and the symmetric (right panel) phase. For comparison we give the suitably 
normalized results from ||5T[ (dashed line). 
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Table 3: Critical exponents from different approaches (A^ = 4) 
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Table 4: Critical amplitude ratios and couplings from different approaches {N = 4). The last three 
entries are non-universal. 
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from the ERG in different approximations also show a larger deviation than was obtained for = 1, 
in the case = 4 our results differ also from the ERG-results by up to 10%. There are a number of 
reasons for this deviation. First, as in the case of = 1 the derivative expansion induces an error in 
the critical exponents. We expect the error from the neglect of the wave-function renormalization 
to be responsible for most of the difference between the values found in the present work and the ones 
given in ||l9[ and |51]. We have also checked in the three-dimensional effective theory and using the local 



polynomial expansion for the effective action how much the choice of a 0-function cutoff in ( p^ affects 
the result (in lowest order in the derivative expansion). In principle, the choice of the cutoff- function 
has no effect on the results if one solves the full RG-equation. Since we have to make approximations, a 
slight cutoff-dependence is expected. In the case we have checked, the results for the critical exponents 
using a sharp cutoff and using an exponential cutoff as usually used in the framework of the ERG differ 
by about 3%. 

There is however an additional contribution in first order in the derivative expansion for A'^ > 1, being 
of the structure Zi^{p)Lp°^ip'^d'^(p^ ip^ . This term should account for much of the remaining difference to 
the exact exponents. 

The same remarks apply to the universal amplitude-ratio Ry- which is compared to the value found 



in |51] in table 4. In view of these uncertainties the values for the exponents and amplitudes found in the 
present work are consistent with the expectations, although they certainly leave room for improvement 
along the lines discussed above. 

Finally, we consider the A^-dependence of the critical temperature to compare with the leading order 
result of the large- A^ expansion. In figure 10 we have plotted our values for g[T = 0) = 0.1 as a function 
of A^, where again the critical temperatures are given in units of ^fpQ (solid line). For comparison, the 
dashed line shows the result that is obtained in leading order in 1/A^ [^] 



Naive perturbation theory yields (|55| ) which is consistent with our results on the level of 1 — 2% for the 
small value of g'T=o = 0.1 choosen here (see figures 5 and 8 for a comparison of with (|55[) as a function 



of the coupling). Also given are the results from the ERG in the imaginary-time formulation [25| as 
squares. Again we observe good agreement as was already seen in figures 5 and 8. It is interesting that 
we can also study the limit A^ ^ straightforwardly. 



4 Conclusions and outlook 

In the present paper we have discussed the finite-temperature phase-transitions of self-interacting scalar 
theories with at vanishing temperature spontaneously broken 0(A^)-symmetry. We have used an imple- 
mentation of the Wilsonian renormalization group-approach for field theory in thermal equilibrium in 
the Schwinger-Keldysh (CTP)-formulation j^. 

We have reviewed and discussed this "thermal renormalization group" (TRG)-approach in section 
2. The approach has a number of advantages as compared to the more common RG-approach in the 
framework of imaginary-time thermal field theory: It allows for the direct computation of real-time 
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Figure 10: Tc in units of ^/po as a function of N for g{T = 0) =0.1 (solid line). For 
comparison results from [^5| (squares), and the result of an expansion in 1/A^ are 
also given. Our results are very well reproduced by the naive perturbative result 
(p5D for all N. 



observables in field theory at finite temperature, avoiding the need for analytical continuation. Thus 
quantities like thermal damping rates and the like are easily accessible and first results may be found in 



the literature |1C]. As the imaginary-time RG approach, the method is non-perturbative and allows for 
the treatment of situations in which perturbation theory fails, which is typically the case near second 
order phase-transitions or in strongly coupled theories. 

Since the real-time formulation of thermal field theory clearly distinguishes thermal and quantum- 
fiuctuations and the TRG treats only the thermal fluctuations, one obtains a direct connection between 
the physical couplings of the theory at vanishing and at non-vanishing temperature and there is no 
need to take into account zero-temperature renormalization scheme dependencies. Furthermore the TRG 
manifestly respects gauge-invariance. 

Finally, after a straightforward generalization of the formalism given in section 2, the TRG may easily 
be applied to field theory at finite density (chemical potential) or to theories with a non-thermal density 
matrix. The method thus offers a nonperturbative approach for the study of a number of interesting 
questions in field theory in a hot and possibly dense environment. 

The main part of the present work is concerned with the application of this formalism to 0{N)- 
symmetric scalar theories in 3 + 1 dimensions at temperatures near to the critical one. We have given 
extensive discussions of the cases = 1 and N = 4:, being especially interesting for a number of reasons. 
The case = 1 in three dimensions corresponds to the well known Ising-model. The electroweak sector 
of the standard model is in the universality-class of the Ising model if the zero-temperature Higgs mass 
is at its critical value ~ 70 GeV f[% . 

We have given the critical temperature in terms of the scalar self-coupling at T = and found good 
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agreement with results found from the exact renormahzation group in the imaginary-time formulation 
(figure 5). The results disagree with the values found in the environmentally friendly RG [29| but 



are consistent with results from the auxiliary mass method |3C], where the latter yields results that are 
basically identical to the perturbative values for couplings up to 5 = 1/3 (figure 3 of the last reference of 
P0[|). Although the approach of [^] differs from the one used in the present work, the discrepancies of 
the results on the critical temperature are somewhat surprising and should be further understood. 

We have then at length discussed the critical behavior of the theory, focusing on the universal aspects 
encoded in critical exponents, amplitude ratios and the scaling equation of state (tables 1 and 2, figures 
6 and 7). Whereas the exponents found in the present work are not as accurate as the results found 
from Borel-resummed e-expansion and other methods working directly in the three-dimensional effective 
model, the scaling equation of state compares very favorably with results from Monte-Carlo simulations 
in three dimensions. This demonstrates the power of the TRG and also demonstrates explicitly in a 
nonperturbative context how the critical behavior of the finite temperature theory is governed by the 
three-dimensional Euclidean model. 

After the discussion of = 1 we have turned to the case A^ = 4, which is presumably in the 
same universality-class as the chiral phase-transition in the case of 2 massless flavors. We find a similar 
situation as in the case of A^ = 1: The critical temperatures are compatible with the results from the 
ERG in imaginary-time ||2^ (figure 8). For this case, results from the environmentally friendly RG or 
the auxiliary mass method are not available. 

Concerning the critical behavior, we again find a scaling form of the equation of state in very good 
agreement with the results from the exact renormahzation group as given in |5l| (figure 9) but a somewhat 
larger difference in the critical exponents and amplitude-ratios (tables 3 and 4). We have argued that 
this situation should improve going to higher order in the derivative-expansion. 

Finally we have given the critical temperature as a function of A^ for a large range of A^, using a small 
value of the zero-temperature coupling gT=o = 0.1 (figure 10). We are able to study the limit A^ ^ 0, 
which is of interest for statistical physics and find results consistent with the values obtained from the 
exact renormahzation group in the imaginary-time formulation and also with naive perturbation theory^. 

Altogether the TRG has proven a useful and flexible tool in the study of field theory at finite temper- 
ature. It allows for the study of universal and non-universal quantities and may be extended to theories 
involving gauge-fields |l^ and fermions or to situations at finite density or more general density- matrices. 
It offers a nonperturbative way of studying dimensional reduction and furthermore allows for the inves- 
tigation of real-time quantities and questions related to the dynamics of systems at finite temperature. 
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